Magnetotransport in the insulating regime of Mn doped GaAs 



00 
O 

o 

(N 



o 

CZ3 



C3 



c 

O 

o 



> 
m 

O 
oo 
O 



X 



Louis-Francois ArsenauuQ 
Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854 

B. MovagharQ P. Desjardins, and A. Yclon 
Departement de Genie Physique and Regroupement Quebecois sur les Materiaux de Pointe (RQMP) 
Ecole Polytechnique de Montreal, CP. 6079, Succursale "Centre-Ville", Montreal (Quebec), H3C 3A7, Canada 

(Dated: July 9, 2008) 

We consider transport in the insulating regime of GaMnAs. We calculate the resistance, mag- 
netoresitance and Hall effect, assuming that the Fermi energy is in the region of localized states 
above the valence band mobility edge. Both hopping and activated band transport contributions 
are included. The anomalous Hall current from band states is very different from the hopping Hall 
current and has extrinsic (skew) and intrinsic (Luttinger) contributions. Comparison with experi- 
ment allows us to assess the degree to which band and hopping contribution determine each of the 
three transport coefficients in a particular temperature range. There are strong indications that 
the insulating state transport in GaMnAs is controlled primarily by extended state, band edge, 
transport rather than by variable range hopping, as reported in the literature. 

PACS numbers: 75.50.Pp, 75.47.-m, 72.20.Ee, 72.15.Rn 



I. INTRODUCTION 

Magnetically doped semiconductors constitute a very 
active field of research, for good reasons. As magnets, 
they are expected to keep some of the useful properties 
of the host system. As doped semiconductors, they have 
their own usual applications now combined with new 
functionality. For example, with Quantum Well (QW) 
multilayers and Quantum Dot (QD) Stranski-Krastanov 
growth technology 1 , it is possible that spins are localized 
in small clusters as in InMnAsS and in QD layers 
of MnGaAsi. This could be used in building small 
magnets, memories and spintronic switches. Therefore, 
localized state transport, magnetism of cluster-localized 
states, and tunneling are very interesting and important 
phenomena. However, we must begin with a clear 
understanding of the bulk materials, whose electronic 
transport properties in the insulating regime we investi- 
gate here. We demonstrate that the observed properties 
cannot be explained solely by hopping in homogeneous 
media. For this purpose, we need to develop the hopping 
theory alongside the band theory and compare their 
predictions. The approach in this paper has a wide 
range of applications and is, we believe, useful to the 
wider magnetism and device communities. The theory 
will enable us to treat resistance, magnetoresistance, 
thermopower, normal and anomalous Hall effect in the 
high resistance regime of DMS. 

Mn doped semiconductors » are disordered alloys, 
so one cannot expect sharp band edges, and indeed, 
one should observe Anderson localized states above 
(for holes) the mobility edges. The long range fer- 
romagnetism is mediated by holes, which lower their 
energy each time they interact with a Manganese spin of 
opposite spin. The mobile valence band p-holes interact 
with the localized moments via the antiferromagnetic 



exchange coupling 



Hpd = 



me{N s } s,s 



J m^ms ) 



(1) 



where J p d is the hole-Mn coupling, cr ss / is a vector 
containing the Pauli's matrices i.e. [<r x ,<t v ,<r z ]. The 
indices, s and s' indicate which terms of the 2x2 matrix 
we are considering. Finally, S m is the Mn spin operator 
at site m and the c ms are creation and annihila- 

tion operators for a carrier of spin s at site m. This 
Hamiltonian, applied to GaMnAs, has been discussed 
in many papers and we refer the reader to the original 
literature (Ref. d and references therein). The transport 
and magnetism was recently subjected to a coherent 
potential approximation (CPA) treatment^. The CPA 
is a mean field theory which can handle delocalized 
states up to the short mean free path or diffusive limit, 
but it does not reproduce localization. Localized states 
and hopping transport have to be treated separately, as 
shown below. 

Near the band edges, the eigenstates of one parti- 
cle systems are localized for energies up to the mobility 
edges The localization starts when the density 

of states is below a critical value determined by the 
diffusivity and localization length^. The nature of the 
mobility edges in disordered magnets is expected to 
be more complex than that of nonmagnetic alloys, but 
many of the usual features should remain valid. In a 
magnet such as GaMnAs, we now have two spin-split 
bands, and two mobility edges. Though conductivity, 
magnetoresistance and Hall coefficient have been mea- 
sured in the insulating limit^, the dependence of the 
mobility edges on magnetism has not been investigated. 
Neither is it clear from the present understanding of the 
data whether all transport coefficients are truly in the 
hopping limit, or whether there is a substantial band 
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edge delocalized contribution which looks like hopping. 

One usually says that transport is in the hopping limit 

when the measured conductivity scales with temperature 
1-1 

as exp — (-y-) 4 ■ Here, Xb is the renormalized Mott 

temperature To = 24 x 2.7 a 3 / (irp(s f)kB) with a denot- 
ing the inverse localization length^, ks the Boltzmann 
constant, and p(sf) is the density of states at the 
Fermi level and 2.7 the percolation factor. But hopping 
transport with Coulomb correlation gaps often scale as 



exp 



T I 



and there are many situations in which 



materials exhibit exp 



'To 
v T 



< n < 1 laws, 



which are obviously not due to variable range hopping 
(VRH) between localized states (see Ref. [^). EmmiS, 
for example, has shown that small polarons will also 
produce such behavior, for reasons not related to VRH. 
Indeed, one can say that almost all disordered insulators 
exhibit such a weaker-than-exponential law. Often it 
is due to granularity and Coulomb correlation o 9 ' 11 . 
Furthermore, it is often not easy to distinguish between 
(i ~ 1/3, 1/4 . . .etc because the data are collected over a 
limited range. So one has to be cautious and examine 
other possible scenarios. This is what we propose to do 
here. 

We do this by presenting a complete theory of 
magneto-transport which is useful for other materials as 
well. But, as an application of the formalism developed 
here, we concentrate our analysis on Mn doped GaAs. 
On the basis of experimental data for magnetoresis- 
tance and thermopower, we question the assumption 
of Mott-hopping in inslating GaMnAs and we offer an 
alternative explanation. Concerning the Hall effect, we 
recall that Allen et alM have shown that the Anomalous 
Hall Effect (AHE) in the insulating regime of GaMnAs 
is not just anomalous in the usual sense because it scales 
with magnetism, but also has an anomalous sign. An 
explanation of this anomaly is apparently contained in 
the hopping Hall effect model of Burkov and Balentai^. 
Here, we propose an alternative theory which includes 
both hopping and band edge conduction. 

This paper is structured as follows: we first present the 
basic definition of the Hall effect in magnets. The pro- 
cesses that explain transport in the insulating limit are 
then discussed. All the main contributions to the Hall 
effect are then introduced, including hopping, intrinsic 
and extrinsic, and the behavior near the mobility edge. 
We then discuss the magnetoresistance and the different 
contributions to it are introduced. Following that, 
all the contributions to the transport propertires are 
summarized and discussed in relation with experiments. 
Finally, we present a discussion of our results. We argue 
that transport in the insulating regime of the magnetic 
semiconductors of the GaMnAs type above 10K, is 
mainly one of delocalized band edge conduction. 



II. HALL EFFECT 

The experimentally measured Hall coefficient Rh is 
often written as 



Rh — Rn 4 
Ra = (apxoi 



Ra 
+ b P 



2 ) 



M z /B e z 



(2) 



where a, b are constants and p xx is the resistivity. B z xt 
is the external magnetic induction, defined as pQB ext , 
where H z xt is the external magnetic field. The first term 
Rn is the normal Hall coefficient and scales with resistiv- 
ity in the usual way and the second Ra is the anomalous 
term which, in general, can have two components, one 
linear and the other quadratic with resistivity^ and is 
proportional to the magnetization M z . The general rela- 
tion for the Hall coefficient is^ 



Rh = 



Pyx 
J^ext 



with 



Pyx — 



-Re{a xy } 



Re {cr X x} Re Wyy} - Re {a yx }Re {a xy } '' 



(3) 



(4) 



where a xy and a xx denote the transverse and normal con- 
ductivity respectively. We may suppose that, as there is 
no voltage applied in the y direction o yx a xy <§; o xx a yy . 
Also, if we consider an isotropic system in the xy plane 
&xx = &yy and thus we may write the Hall coefficient 



Rh 



Re{a xy (B z )} 
B% xt [Re {a xx {B z )}\ 



(5) 



Karplus and Luttinger 15 pointed out that the B-field in- 
volves the magnetic moment of the material via the in- 
ternal magnetization M z 

B z = po [H e z xt + (1 - N)M Z ] ee B cxt + /Xo(l - N)M„ 

(6) 

where N is the demagnetizing factor. Thus, the magne- 
tization term is implicit in the normal contribution as a 
shift in the magnetic field. However this form is not nor- 
mally sufficient to explain the much larger magnetization 
contribution observed in ferromagnets 15 . Karplus and 
Luttinger then developed the first theory of the intrin- 
sic AHEi£. Throughout, we shall, for simplicity, suppose 
a thin film with the magnetic field perpendicular to the 
plane and thus, we take N — 1, unless otherwise men- 
tioned. Let us now consider the problem of conduction. 



III. HALL EFFECT AND CONDUCTIVITY IN 
THE INSULATING LIMIT 

When the Fermi level is at the band edge, Anderson 
localization sets in, and Fermi level transport cannot 
be described using only extended states. The very low 
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temperature transport process is by hopping, as shown 
by Allen et al& for GaMnAs. These authors have also 
measured the Hall effect and, as already mentioned 
above, observed that there is a sign anomaly for the 
AHE in this regime. They argue that the transport 
is by hopping from localized level to localized level 
and not by thermally assisted hole ionization into the 
valence band. The magnetoresistance in the insulating 
regime has been measured by Van Esch et air . These 
authors constructed a model to explain the very large 
magnetoresistance they observed, which we will discuss 
later. But no one seems to have yet made a model which 
consistently describes resistance, magnetoresistance, and 
Hall effect, in the insulating regime of GaMnAs. 

The ordinary phonon assisted hopping Hall effect 
has already been treated rigorously and in great 
detai l 16 ! 17 ' 18 . The extraordinary contribution reported 
recently by Allen et al& in Mn doped GaAs, can be 
modeled as follows. We consider p-holes in the valence 
band, spin-orbit coupled as in the models of Baldareshi 
and Liparii^ and Fiete et al££. The holes are bound to 
a negatively charged, spin 5/2, Mn ion. Fiete et aZ/22 
also include the AF exchange coupling of the hole spin 
to the d-spin. 

The Hamiltonian describing the dynamics of holes 
localized around the Mn sites and which incorporates 
the spin-orbit coupling to first order is given by 



Z,3 



Xi(i,s\li ■ <r\i, s')q s Cj. 



(7) 



- £ Ms (k,z + gci.z) (B z + h z 



In Eq. ([7]), the summation indices include the band and 
site indices, i.e. j = j, 7 where 7 denotes the orbital, 
and we again assume N = 1. The Ei S are diagonal site 
energies. Now consider, the second term, where Uj(B z + 
B so ) is the site-to-site transfer term. The magnetic field 
is incorporated in the Peierls phase. We recall that in the 
presence of the magnetic field we have the usual Peierls 
phase factor 



t n 



B cxt -(R„xR„ 



(8) 



The non-local spin-orbit terms are not specifically 
included in Eq. ((?]) but can be reduced, in first order, 
to a shift of the total magnetic field by an amount B so , 
that is included in tij and which we consider in detail 
below. The atomic spin-orbit coupling Ai, in the third 
term in Eq. l[7]), mixes angular momentum 1 and spin 
bands on the same atom. Assuming, for the sake of 
argument, that we have orbital states which are p-states 
with l z — {1,0,-1} or combinations thereof, then the 
coupling will flip, for example, 1 1,-1/2) — > | — 1,1/2) and 



back again. The z-component of the atomic spin-orbit 
term ( — AZj z a z ) can be included in the sum of the 
Zeeman and Curie Weiss spin splitting energy for the 
carriers, where \xb is the Bohr magneton. This sum is 
denoted by h ZtSp i n in the fourth term of Eq. ([7]). 

The spin-orbit Hamiltonian in tight-binding has 
been examined by Pareek and Brunch for disorder 
induced coupling, by Movaghar and Cochrane 1 ^ for lon- 
gitudinal and Hall transport and in CPA by Arsenault 
et al£*£. In weak spin-orbit coupling, the effect can 
be incorporated into a two-site overlap phase similar to 
the Peierls phase given by Eq. ([8| 4 ' 14 ' 22 . In the hopping 
limit, with no external field, the electric field at a site 
i (multibands are implicit if necessary) is given by the 
sum of the fields emanating from all neighbors 



eZ n (r - R„)(l - p n ) 



47T££o|r — R r , 



(9) 



In Eq. ^j, e is the dielectric constant which screens the 
acceptor site potential, Z n is the effective charge, and the 
factor (1 — p n ) is the probability that the site n is empty 
and therefore charged. In the i to j jump, the third closest 
site in the pathway will contribute an electric field which 
produces a net magnetic field. So the sum in Eq. ([9]) 
can be approximated to include the best neighbor to i 
(see Refs. @, [13, and refs therein, for the definition of 
best neighbor). We can now include this spin-orbit phase 
in the magnetic phase by adding to B ext in Eq. ([8]) the 
simplified effective spin-orbit field 



B so (i -> j) 



2mc 2 



£ 



Z„(l -p n ) 



47T££ |r - RJ 



where the nearest neighbor sum excludes the start site 
i and end site j. A more complete derivation for the 
effective spin-orbit field is presented, for the case with no 
screening (no e in the equation) and where all potentials 
contribute (no 1 — p n in the equation), in Ref. @. Note 
that the integral (i\ | r _^ p |z)is strictly speaking not 
convergent. But, in practice, the orbit radius is never al- 
lowed to be smaller than the effective atomic orbit of the 
valence state so that the cubic singularity does not occur. 

Before determining the various contributions to the 
Hall conductivity, let us examine the effective masses 
which may appear in the expressions which follow, so 
that the meaning is clear. There are three effective 
masses in the problem: 

1. the Kane-Luttinger mass (m*) which controls the 
kinetic energy in the starting Hamiltonian (effective 
mass Hamiltonian); 

2. the spin-orbit mass which is equal to the bare elec- 
tron mass 
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3. an effective mass which is a result of the sum rule, 
as discussed in Refs. andU^, and which can be- 
come the diffusivity in strong disorder but is the 
Kane-Luttinger mass in weak disorder. 



Here p xy is the Hall resistivity, ao the conductivity 
prefactor and we have used N = 1. Note that for the 
remainder of the paper we will no longer distinguish the 
real part of the conductivity from total conductivity, 
since they are the same at zero frequency. 



A. Hopping contribution 

The normal Hall term arising from the external B = 
B z term has been derived by the hopping interference 
method in the pure diffusive limit, and in the VRH 
limi t 12 ' 17 ' 26 and applied to the case of phonon assisted 
hopping and percolation in a one band model. Thus, 
we may add the spin-orbit magnetic field to the external 
field and carry out the configurational average as before. 
So we replace B z by 



Bl xt + fx (l - N)M Z + B so , 



where 



B s 



E 



Z n (l -Pn) 



2mc 2 *—f 47T££ol r — R-nl 3 




(11) 



(12) 



where (a z ) is the polarization of the holes and M z the 
magnetization of the sample. The analysis can now be 
done as previously. The high density hopping system 
for example, has been treated analytically in Ref. 
and the phonon assisted percolation regime in Refs. 
El and m 



The second term in Eq. (jlip depends upon the geometry 
of the sample and is zero in thin films under normal 
fields^. The sum in the third term runs over the third 
neighbor in the triad and is configurationally averaged to 
the optimum second nearest neighbor distance (as shown 

in Ref. GJ), which is |R ( - R„| = (H a (T)) cx (2a.)*. 
Here, To is the renormalized Mott temperature, defined 
previously. 

Replacing Eq. (Ti"Tj) in the expression of Ref. al- 
lows us to calculate the Hall mobility (sign is negative 
for electrons and positive for holes) and conductivity 
as was done before^. Assuming a constant density of 
states, we obtain 



Pxy 



ReWxy} _^_§(l«)i 



5f'Re{ff ra } 
Re{<T xx } = foe" 



and 



Pxy 

Bf* ~ Re{a xx } 



Br 



l 



hop 



Re{pxx} ' 



B' 



-A 



hop 



(13) 
(14) 

(15) 



The factor 3/8 in Eq. (TT3"]) is a consequence of averaging 

i 

over the third site in the triad (2aR 3 ) ~ M(lb.) 4 . 
When the density of states is quadratic in energy, rather 

than constant, we replace (-^) 4 by (^#) 2 . This is also 
true for the Coulomb gap hopping proble m 16 ! 23 ' 27 ! 28 . 



" ) 



Allen et alfi- observed an experimental exp 
dependence of conductivity. This is characteristic of 
a Coulomb gap or hopping in a quadratic density of 

states. They also found a near p xx behavior of the Hall 
resistivity. This is also reproduced by the present theory 
as shown in Eq. (|15[) . Finally, they found normal and 
anomalous processes which have opposite signs. The 
sign anomaly can be explained by the sign of {a z } which 
can be opposite to the overall magnetization since it 
refers only to carriers at the Fermi level. 

The maximum strength of the AHE predicted by 
Eqs. (TT3"]) . (TT4"]) and (fTS"]) is determined by the magni- 
tude Xhop- An estimate which includes the factor 1000 
reduction due to the cube of the average distance to the 
3 rd neighbor, gives ~ 3 x 10 3 Gauss. This spin-orbit 
term is somewhat too small at T = 25K when the result 
is compared to the experimental data of Ref. H|. The 
data suggest that the AHE is more important than the 
Normal Hall effect, even at an external field of 10T. 
The measured order of magnitude is close to the theory 
estimate only in the truly very low temperature hopping 
regime where AHE ~ NHE i.e when T ~ 10K or less. 
The ratio of AHE to NHE is always roughly 5 or 6 when 
we look at the data in the metallic regime. It seems 
therefore that though the simple hopping model may 
be right for the truly low temperature hopping regime 
T < 10K, at higher temperatures the spin-orbit electrical 
fields which are needed to produce the sideways force 
are larger than those predicted by hopping theory. 

This is possibly because, in this regime, the Hall 
charge transport is, in fact, already activated to above 
the mobility edge and is not really hopping-like. We shall 
examine this possibility in the next section. This view is 
also supported by magnetoresistance 7 - and thermopower 
data^(although Osinniy et aL— say otherwise). Thus, 
the picture of the third site in the triad providing the 
spin-orbit field is attractive, but the effect produced 
by this mechanism is probably not large enough to 
explain the quantitative values. States excited above the 
mobility edge undergo skew scattering and intrinsic Hall 
scattering and feel the full effect of the lattice-induced 
spin-orbit enhancement, as discussed by Chazalviel 3 -^, 
Berger- 3 ™ Fivaz^S, and Engel et al. 33 . The previous work 
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of De Andrada et alM', for example, gives a maximum 
enhancement factor of the bare coupling of the form 



"fenh 



(4) 



A(2E g 



E g {E g + 3A) (3E g + A) 



(16) 



where A is the Kane spin-orbit coupling energy, E g is the 
energy gap of the semiconductor and m* is the effective 
mass. One can see that in GaAs, the enhancement of the 
spin-orbit coupling is considerable. Thus, the pathways 
above the mobility edge, even if not necessarily always 
dominant for the ordinary transport, can be expected to 
be dominant for the AHE. It is difficult to see where such 
an enhancement would come from in hopping conduction. 
Including the many orbital bands will not strongly affect 
the results in the hopping limit. The marginal but im- 
portant contributions from the hopping pathways forces 
us to now look for other possible explanations, or other 
contributions to the hopping channel specially at higher 
temperatures. 



In the preceding equations, hb = f^, fi^a) is the Fermi 
function, n as is the number of electrons in the state a 
with spin s and Ag^ z is the zz component of the (/-shift 
tensor. This component can be related to the result first 
derived by Karplus and Luttinger— by writing approxi- 
mately Ag ~ X/ Agap where A sap is a (Bloch) subband 
gap and A the spin-orbit coupling strength. The g-shift 
is used here as in the language of Chazalviel^, and in- 
cludes the Bloch band enhancement. Note that the in- 
trinsic contribution given by Eq. (TTTj) apparently docs 
not scale with the conductivity. But this is only true in a 
highly ordered medium in which the scattering lincwidth 
is smaller than the intersubband energies. In this limit, 
it only appears as a lifetime term in a basis of well de- 
fined Bloch states. In a strongly disordered system, there 
are no well defined Kohn-Luttinger bands as such, and 
the g-shift will contain the relaxation time term in the 
denominator as well. When the linewidth is larger than 
the dominant excitation energies in the sum, then the 
relaxation time will enter the g-shift and the Luttingcr 
term will also scale with resistance. 



B. Intrinsic Luttinger contributions from the 
region above the mobility edge 

Consider the contributions to transport coming from 
excited states below (for holes) the mobility edge. For 
delocalized states, there are skew scattering Hall contri- 
butions due to impurity potentials. This has been treated 
in Refs. [33| and l35l . There is also an intrinsic Luttinger 
contribution which can be written, in the language of 
Chazalviel 3 -^, and as rederived to first order in spin orbit 
coupling by Arsenault and Movaghar— . We write the 
transverse conductivity as 



■xy 



e 

TT e 



df(e a , 
de a 



(17) 



where 



C. Magnetic field dependent resistance and Hall 
effect contributions from the region above the 
mobility edge 



n a \ - n a i 



(18) 



e = — |e| for electrons and e = \e\ for holes. Finally, the 
zz component of the effective (/-shift tensor is 



A<7? = 

E 



4m 2 c 2 



(W(r) x p) 



e/3 



(P\L z \a) 



We can now give a more complete formula which in- 
cludes the hopping channel and the mobility edge con- 
tributions and which allows us to predict the magnetore- 
sistance as well. From standard mobility edge transport 
theory£i2&, and including the normal Hall effect, the skew 



scattering contributio n 33 ^ 5 
Luttinger term, we have 



and the intrinsic Karplus- 



'xy 



de (-^f) P{c)D dlff {e) (^k^) + + ^vM P ^ g (T) + 4 y , 



(20) 



where the so called intrinsic AHE is 



and 



,/:|-^]p( E )AA 5 (^(e) (21) ^x = e 2 



dfhje) 
de 



p(e)D diff (e) + a hopp i n g (T) . 

(22) 
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The second part of the first term in Eq. (|20|) is the skew 
scattering contribution with — denoting the skew scat- 
tering rate. Here, //i(e) is the hole distribution. In 
Eq. IJ3DJ and Eq. Qg), D dtff (e) is the band diffusivity 
at energy e, and by localization theory 3 ^, is of the form 



Ddiff{e) = D 



(23) 



near the mobility edge. Note that the Luttinger term 
contains the quantum diffusivity h/m* and not the ac- 
tual diffusivity Ddiff. In Eq. (T23|) D ~ lcm 2 /s, is 
the density of states for which one normally adopts an 
exponential form: 



p(e) = p e 



(24) 



with < — e < —£ c , where po is the value at the Fermi 
level and Eq measures the exponential steepness. The 
spin-orbit coupling is contained in the g-shift factor. If 
the intersubband gaps A gap is definable, the (/-shift is 
Ag a ~ A A as in Karplus and Luttinger 15 . In highly 
disordered alloys, there are no well defined subband gaps, 
as assumed by Karplus and Luttinger^ 5 . and Jungwirth 
et al. 3 . Then, one has to evaluate the g-shift more 
rigorously. The g-value of the Mn spin + hole complex 
(~1.98) has been measured by ferromagnetic resonance 
and is in the range (-0.05 to -0.2), hole like (Ag < 0), and 
increases with Mn-hole doping (see the review by Liu 
and Furdyna 37 ). Note that the sign anomaly can again, 
as in hopping, be due to the fact that the holes at the 
Fermi level are polarized opposite to the bulk magnetism. 

The hole Fermi function //t(e) is, to a good approxima- 
tion, a Boltzmann function when \et — e c \ ksT. The 
mobility edge for holes is denoted by e c , and the energy 
dependence of p(e)Ddif /(e) dominates the longitudinal 
conductivity. The anomalous skew scattering Hall 
contribution will be dominated by the spin-orbit band 
term via in Eq. (|20"|) 33 - 35 . The skew scattering 

contribution is not negligible because of the very large 
spin-orbit coupling enhancement (10 4 ), implicit in 
Eq. (fi"6|) , when applied to GaMnAs. It is useful to recall 
that Engel et ali~ have estimated the skew factor to be 
t" ~ 9M f° r doped GaAs. Thus, when comparing 
-j-(o~ z ) with uj c t, and noting that r ~ jf^zjj, e <e c (for 
holes) , we see that the skew scattering band contribution 
to the AHE is comparable to the ordinary Hall effect 
even near band edges. The Luttinger term (Eq. (J2TJ) ) 
is ~ {10~ 4 — 10~ 3 } x Ag(a z ) and will dominate when 
Ag > 10 -3 . We return to this later. The mobility 
edge channel gives rise to very concrete formulae for the 
behavior of the magnetoresistance and Hall effect. So 
now we should examine the hopping magnetoresistance 
and then we may be able to decide which mechanism 
dominates in each range of B-field and temperature. 



IV. MAGNETORESISTANCE 
A. Hopping Magnetoresistance 

The relative importance of the hopping and band edge 
weak localization contributions could, in principle, be 
tested by measuring the corresponding magnetoresis- 
tance. We shall now show that these are very different, 
both in structure and in order of magnitude. It is surpris- 
ing that magnetoresistance was not measured by Allen 
et o/4. The magnetoresitance in the apparently hopping 
like regime was measured in a similarly low doped (2%) 
system by Van Esch et alJ-, who found a strong negative 
magnetoresistance, which they interpreted as being due 
to the orbital expansion of the localized states in a 
B-field. This, the authors suggested, is caused by the 
competition between the spin Zeeman energy and the 
antiferromagnetic hole-Mn coupling. The external field 
is, however, far too small to interfere with the strong an- 
tiferromagnetic J p d coupling of ~ 0.5eV— . The Van Esch 
et alX magnetoresistance data actually look very similar 
to those obtained for metallic films by Matsukura et al£^. 

Let us now examine the complete hopping magne- 
toresistance and then evaluate whether it can explain 
the available data£ in the insulating regime. The theory 
will also allow us to make predictions concerning what 
should be the behavior of the hopping magnetoresistance. 

Specifically, hopping spin and orbital magnetoresis- 
tance (MR) do exis t 39 ' 40 i 41 , and have been discussed by 
Movaghar and Schweitzer— as applied to a-Si. Hopping 
magnetoresistance usually has a positive contribution 
from orbital shrinking of wavefunctions at very high 
field, B > IT (Mikoshiba effect). This is a one-body 
effect. Localized states also have a finite Hubbard U, 
so there is also a two-body contribution to the spin 
hopping magnetoresistance. An electron from below 
the Fermi level can, in general, hop to an empty state, 
or one which is already occupied by another electron 
below. The doubly occupied level is shifted up by the 
Hubbard U and, if U is not too large, this takes some 
of the doubly occupied states to the Fermi level. The 
jump to the occupied state takes place in the triplet or 
singlet configuration. The triplet hop being forbidden, 
the three triplet states would have to first convert into 
singlets. This process depends sensitively on (small) 
transverse (B x ,B y ) magnetic fields within the system^ 2 -. 
The triplet-to-singlet transformation can also occur 
during the jump via the spin-orbit coupling Al • cr. We 
can write the totality of effects as 4 ^ 



r h x °/(B z ,T) = a e'^y 4 



N e + N s + N T 



.(■ 



(WcT) 2 + 1 



(25) 
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In Eq. (j25|) . £ is a quantity that will be given below. 
N e , Ns and are probabilities that the hop is to a 
previously empty, singlet or triplet state. Ns and Nt 
depend upon the magnetic field as N s — and 

Nt = 3 4 . The direct spin-orbit assisted triplet-to- 
singlet hop is proportional to ^ where v Q is the hop- 
ping (singlet) prefactor. The internal random fluctuat- 
ing magnetic fields Bi nt (t) which allow spin admixtures 
to occur are mainly due to fluctuations in g xz , g yz and 
the anisotropic part of the dipole-dipole interaction. For 
weak coupling, we have 



( B int) = ( B dip) + 



l r 



&9yzB e z 



(26) 



The dipolar internal fields {B 2 ^) 1 ^ 2 are of the order of 
the ESR linewidths smaller than lOmT, the g-admixture 
is ~ 1CP 3 . Saturation occurs when the frequencies 
corresponding to Eq. (|26p are comparable with typical 
hopping frequencies. A similar model has been used 
recently to explain the magncto-elcctro-lumincscence in 
polymers^. 

The Mikoshiba effect, mentioned earlier, is spin in- 
dependent and changes the phonon assisted hopping 
rate as 



Wi„ 



-2aR ir 



(27) 



averaged over the angle between the B-ficld and R 
and with Sj > £j. At optimum hopping length, i?f 7 

goes to ((aR) 3 } = (f)^^) 1 and 



,1 



goes to 



I (lb} 
4 V T ) 



Inserting these terms in Eq. (|27p gives us 

L 3 

i _J_ in Eq pgj) t lowest order, we 



(To\2 (3\4 1 

\ T J Is/ to 1 

may neglect the change in T caused by B z . There is 
also an increase of the energy of the localized state. This 
shift lowers the ionization energy to the mobility edge 
which, in first order perturbation theory is 



e l {B z 



(eB z ) 2 
2a 2 m 



(28) 



where a, is the inverse localization length of the i th level. 
Finally, there is also a small quantum interference magne- 
toresistance contribution due to the fact that the hopping 
carrier could also virtually visit the second best neighbor 
site on the way to its best neighbor. This contribution 
has been analyzed in detail by Schirmacher 44 . The spa- 
tial part of the transfer rate is, to first order, proportional 
to 



\T iS (B) 



+0 e -Tf B ( R 'X R i) 



E 

i 



il Ij 
— ( 

£l 



fB-(RiXR|+R[XRj) 



(29) 



This now includes the direct transfer, and the inter- 
ference term from a third site. This means that the 



Miller- Abrahams 4 ^ hopping spatial factor, \t%\ 2 , should 
be replaced by Eq. (|2"9")) which is a small but magnetic 
field dependent factor. One need take only one term in 
the I sum, the second nearest neighbor term. The opti- 
mized second neighbor 2aRu and El k ~ B '^ wm a ^ so sca l e as 

7 (S-) 4 but with a constant factor 7, which instead of 
the first neighbor values, 3/4 and 1/4 respectively, will 
be somewhat larger, ~ 39/32 for the spatial length and 
13/32 for the energy. We have already used this for the 
hopping Hall effect (Eq. (flU)) ). The sign depends upon 
the geometry of the triad and can be positive or negative. 



We conclude that in 

2 



a by cr 



Eq. (f25|) , we have to replace 
where () is the spatial/energy 



average over the best second neighbor (see Refs. [17j, 
I23I and Refs therein for the definition of best second 
neighbor). Keeping the best third neighbor I and 
assuming Ri = 0, we obtain a quantum interference 
contribution 



Aa (B) = cr 



Uj{£i — Si) 



COS 



(|B • (R ; x R,)) - 1 



(30) 

The average in Eq. (|30p is over the percolation network 
and gives a negative magnetoresistance, according to 
Schirmacher 44 . In general, one has to note that all hop- 
ping contributions below IT are at best a few per cent, 
and nowhere near the large values observed by Van Esch 
et alZ. This completes the hopping analysis. Now con- 
sider the contribution to the magnetoresistance coming 
from excitations to the band edge. 



B. The orbital magnetoresistance from the region 
near the mobility edge 

Following Khmelitskii et al^ and Movaghar and 
Roth 46 , we have for holes and small B-fields. the follow- 
ing shift for the mobility edge: 



\e f (B z )- Sc (B z )\~ £f (0)-e c (0) 



1 



hiu c 



MO), 

(31) 

where k is a constant of order 1 and tu c is the cyclotron 
frequency of the holes. In Eq. (|3Tj) . we have neglected 
the magnetic field shift in the Fermi energy compared 
to that of the mobility edge. This result is valid for 
fields at which ftw c < h/r, cyclotron energy is less than 
the scattering broadening, and gives, from Eq. pip , an 



exp 



n[hu c e c (0)] 2 
k B T 



increase in the conductivity over 



limited field range. 

The above describes the well known orbital der- 
ealization effect. The fully self-consistent shift in the 
effective activation energy, i.e. the exact value of the 
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one-body £f{B z ) — e c (B z ) shift, which includes the 
ferromagnetic effects, and high B-fields is not known 
and this would need a detailed study of localization in a 
disordered ferromagnet. However, it can be inferred that 
in the present case, field induced magnetic alignment 
is going to lower disorder and thus increase the range 
of delocalized states by shifting e c (B z , M z ) up with M z 
and Ef(B z ) down with B z . At low temperatures, the 
magnetism can be saturated at relatively low fields and 
will not change the electronic structure very much below 



V. SUMMARY OF CONTRIBUTIONS, 
INCLUDING SPIN DISORDER SCATTERING 
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We now have a theory which allows us to calculate, 
in principle, the resistance, magnetoresistance, and 
Hall coefficient when the Fermi level is in the region of 
localized states. We have also included the band region 
at, or just above, the mobility edge including the spin 
scattering effects which can be described by the CPA 
(Coherent Potential Approximation). This is discussed 
in Appendix [C] The Hall conductivity has already been 
given by Eq. |20|) . 

Consider first the complete longitudinal conductiv- 
ity in the presence of a magnetic field B z and including 
orbital and spin contributions from localized and 
delocalized states, this is given by 



7T 



de 



Sb 



de 



Pedge{s)Do 



e - e c (B z ) 



e c {B z ) 



d£ I -7^- ) Pband(£)D C PA(£,B z ) 

hop i 



(32) 



where ej, is the energy above which the CPA stops being 
applicable. This energy is close to, but a little bit below, 
the mobility edge (we are considering hole) . This expres- 
sion includes the band edge (first term), the high mobility 
band term (second term) , and the third term, which is the 
hopping contribution. In the first, or band edge term, the 
magnetic field dependence is basically due to the B-field 
induced derealization, and not to spin. In the second, 
or band term, which takes over at higher temperatures, 
the magnetoresistance is mainly due to the spin-disorder 
scattering, and is negative, because the magnetic field re- 
duces the disorder, and thus reduces the scattering rate. 
Figs. Q] and [5] are examples of this behavior. Unfortu- 
nately, we do not have the analytical form for this term, 
but we have the numerical CPA results shown in Fig. [1] 
and Fig. O Also, from Eq. (15) of Ref. 0, we can create 
an approximation to the CPA from perturbation theory. 



FIG. 1: Relative magnetoresistivity in the CPA for carrier 
concentration p = O.lx for impurity concentration x — 0.053 
and J p d — 0.35W , with T as a parameter. W correspond to 
half the bandwidth. Figure taken from Ref. A. This is the 
expected behavior just below (hole) the mobility edge 




0.005 0.01 
kgT/W 



0.015 



FIG. 2: Resistivity as a function of temperature for various 
magnetic fields. Others parameters are x = 0.053, Enm = 
Em = 0, Jpd — 0.35W and p = 0.1a;. W correspond to half 
the bandwidth. The solid line is the zero field result, the 
dashed line is for ^lbB = 1x10" W, the dotted line is for 
HbB = lxl0 -4 W and the dashed-dotted line is for (ibB = 
3xl0" 4 iy . Figure taken from Ref. QJ 



By defining, the pd scattering time as 



1 _ JpdPband(e) 

T p d h 



.2 (M t 



,,;^v . ( ZT C M Z 

_ M _j tanh ^____ 



(33) 
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the CPA diffusivity can be written 



D CPA ~ v(ef 



(34) 



where r e is the non magnetic relaxation time and 

v{e) — J J^r is the band velocity at energy e. Eqs. (f3"3")l 

and (f34|) could be interpreted as a first semi-analytical 
approximation to the CPA result (see Fig. [1] for a 
numerical evaluation). The spin-hopping term is given 
by (J25|). 

The density of states has been divided into three 
sections corresponding to the region near the Fermi level 
(hopping), the region at the band edge (Eq. (HH)), and 
the delocalized region just below e c , in the CPA regime. 

We have established that when the Fermi level is 
in the region of localized states, both band edge and 



hopping will, in general, contribute to transport. The 
question is to what extend is it one or the other? This 
obviously depends upon temperature. But it has been 
argued that if the observed conductivity scales as Mott's 
law, then the explanation cannot be anything else but 
VRH (variable range hopping), for all the transport 
coefficients. There are some facts that speak against 
this simple interpretation. The thermopower data are 
not hopping-like^. Our estimates show that the AHE 8 , 
can, at very best, only be due to hopping at very low 
temperatures. But the Hall data, thermopower and 
magnetoresistance are all on different samples, so one 
has to be very careful and rigorous, and examine all 
potential contributions. 

Indeed, with the usual assumption of exponential 
density of states (Eq. ([24])), the first term in Eq. (f32|) 
will give an activated band edge transport behavior, so, 
in this limit, we have the complete result, 



J 



a xx = e 2 D Q p Q (ef(B z )) 



JJ „ fe. „ T ""I" ffn 



e c (B z ) 
*Z x E (T,B z ) + a h x °/(T,B z 



^° X P (T,B Z ) 



(35) 



Using Eq. (|21|) and replacing the quantum diffusivity 
h/m* with the normal diffusivity, we have 



J 



BE 



(T,B Z ) 



eB zT {s c {B z )) 



Ag(e c (B z ))(a z ) £ 
I 



° h ° P {T,B z ). 



(36) 



If we had kept the original form of Eq. (|21j) , the depen- 
dence of the Hall coefficient (Rh) on the resistance would 
be resistance squared as observed in most ferromagnets. 



To complete the analysis we also need Eq. (|3Tj) . The 
conclusion one can draw by looking at Eqs. (l35|) . ([36]) and 
(f31~j) is that, as temperature increases, an activated law 

should be observed, with a large exp (k£ c -\/ fiu) c /e c /kBT\ 

increase in conductance with B (negative magnetoresis- 
tance). 



A word of caution: the shift in mobility edge with 
B-field as applied to electrons or holes assumes that the 
excited state wavefunction has time to build an Ander- 
son level with all its features^. This is not necessarily 
true for highly excited states. As an example, take 
a-Si:H where no such negative MR is observed despite 



activation above the mobility edge. In the present case, 
the mobility edge is very close to the Fermi level so that 
it is reasonable to assume that this mechanism is also 
operating. 

/ \ 1/3 

Van Esch et al. 7 fit their data to a ( J law. 

But most often a Mott law is observed in conduction. 
Osinniy et al£^ have analyzed the thermopower in low 
doped insulating samples and observed an activated law 
(see Appendix IBl for thermopower). The resistance they 
measured also looks activated, although the authors 
do not explicitly state it. So their data would agree 
with the band edge formulae. From this, one would 
have to conclude that both band edge transport and 
hopping is observed, depending upon the sample and 
the transport coefficient. But when we examine the 
order of magnitude of the hopping magnetoresistance we 
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have, at most, a few percent, even at very high fields. 
The data of Van Esch et alZ on magnetoresistance give 
us a strong clue as to what might be happening. They 
have a few hundred percent MR but, with Mott like 
conductivity. So, finally, we have one last question to 
examine. We need to investigate what other mechanism 
can give us a Mott-like law for conduction but a giant 
magnetic band- like magnetoresistance and AHE. Let 
us examine this possibility before we come to a final 
decision on the mechanism. 



A. Mott law from correlation holes or granularity 

The influence of correlations was treated by Anglada 
et a/i^, and more recently, by Dunford et alr^, in the 
context of transport in polymers, granular systems and 
nanoparticle composites. The point is that in highly cor- 
related systems, the density of excited states at the band 
edge need not, and in general will not, be exponential. 
Coulomb energies and spin fluctuations can modify the 
excited state spectrum. Though the exponential struc- 
ture for p(e) is reasonably well established in systems 
which are weakly correlated, such as nonmagnetic amor- 
phous semiconductors 4 ^, one knows that in electron and 
spin glasses, this is not the case, and final state inter- 
actions are important. When an empty state just above 
the Fermi level is filled by an electron excited from the 
valence band, the delocalized hole created in the valence 
band is not immediately neutralized. The hole which is 



created recombines elsewhere in the material, with an 
electron from below the Fermi level, and in this way, 
charge has been transported. The process of the elec- 
tron separating itself from its hole is complex in a highly 
correlated spin system such as GaMnAs. Here, it is the 
carriers which cause the spin alignment in the first place. 
We can expect the density of states for excitations in 
such a system to be more sharply cut off at small energies. 
We shall see that the proposed structure will reproduce a 
coulomb glass type behavior with temperature 4 ^. A simi- 
lar situation is encountered in granular systems^dS, where 
Mott-like laws are observed as a result of inter-grain ther- 
mally or fluctuation-assisted tunneling 49 . Granularity is 
also expected in these materials. Indeed we can simulate 
this behavior with a density of states per volume of the 
form 



Pedge{£) = PbandS W-'+'l-' , (37) 

characteristic of correlated or granular systems, with 
< v < 1 and where e g , si, are constants which measure 
steepness and the density of states at the Fermi level. 
Finally, pband is a constant. 

Knowing the density of states at the band edge, 
we can now proceed to evaluate the Hall effect, longitu- 
dinal conductivity and magnetoresistance at the band 
edge. To see how we can get a Mott type law from 
extended states, we use the band edge term of Eq. (|32|) 
and Eq. (j3Zj) to get 



depbandDo 



\*f<.B*)-'\ 



-(• 



(38) 



38]) at low T gives 
law (see Appendix lAf . 



A steepest descent evaluation of Eq 

for v — 1/2, an exp — (-^) 

Here, the Mott-like laws are a result of the pseudo-gap 
in the density of states or a depletion of low energy exci- 
tations created by long range correlations. 



band edge conductivity with a correlation hole density 
of states evaluated using steepest descent (Appendix [X| , 
and the usual hopping conductivity, would give us (with 
i/ = l/3) 



B. Quantitative comparison: VRH versus 
delocalized band edge conduction 

To asses the relative importance of the various chan- 
nels, let us compare some values starting with a xx . The 



.(band) 



e c {B z )k B T 



DoPba 



[e f (B z )~e c (B z )\ 



3 2\k B T J 



(39) 
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where 



With Eq. (|3ip . for holes, we obtain 



£max = {vk B Te v ) 



(40) 



o~ xx (band) 



c(0)fc s T 



D Pb 



and 



e/(0)-e c (0) 1 + k 



(41) 



and 



Acr X2 



25 



rnax^c 



[e f -e c (0)})' 



(42) 



For the Hall conductivity, we obtain, to a good approxi- 
mation (y = 1/3), for holes 



a xy (band) 
o~ xx (band) 



+ (Ag(e nm.r 



(43) 



Whereas for VRH hopping (without the magnetoresis- 
tance terms) we have 



a xx (hopping) = -^P( £ f) 



3 (To 

4 I T 



1/4- 



-(^) 1/4 



(44) 

It is understood that the T 1 / 4 factors are chosen to 
be the same so that this determines the density of 
states parameters e g . Thus, (3/2) 4 e g /kB — T . The 
quantities, p band ~ 10 20 - 10 21 /cm 3 eV, D ~ 10" 4 m 2 /s, 
p(s f ) ~ 10 17 - 10 19 /cm 3 eV. T or e g is fixed by ex- 
periment. From the data of Allen et al&, we extract 
roughly T = 2.4 x 10 5 , T = 24 x 2.7a 3 /k B p(e /). The 
localization radius corresponds to a binding energy of 
110 meV which is ~ 3nra for a light hole. The corre- 
sponding p{ef) ~ 3 x 10 19 /cm 3 eV extracted is a very 
large number and normally above the value required for 
derealization. 

One can see that the band pseudo-hopping expres- 
sion can easily account for the dc conduction. We 

see that by selecting v, one can create an e v T > 

-(In\ 1/4 

or e \ T I behavior without variable range hopping. 
Thus, if some experiments give one or the other law, 
it only means that the effective density of states is 
different. This has nothing to do with VRH. Further- 
more, from Eqs. (|31[) and (|4"2")l . it gives the correct 
order of magnitude, the right sign and nearly the exact 
form of the magnetoresistance, whereas the hopping 
magnetoresistance is nowhere near the right order of 
magnitude and structure. The magnetism and large 



magnetoresistance observed by Van Esch et alZ at 
very low temperature suggest that inhomogeneity and 
clustering play an important role. Hopping magne- 
toresistance would, however, seem to be the correct 
explanation for the magneto-electroluminescence data 
in the work of Mermer et al£Q where effects are a few 
percent for B-fields of mT, as in the magnetoresistance 
and photoconduction of a:S i 39 ' 42 i 51 . 

Finally, the same optimization applied to the Hall 
effect gives, from Eq. (|4"3")l 



Rh 



max 



B z 



1 



r xx (band) ' 

In this form, we can see that the anomalous term eas- 
ily competes with the normal term since t ~ 10~ 14 s, 
(Ag) - 10~ 2 - lO" 1 and (cr 2 ) £mox - 0.1 - 1. Such a large 
effect is very difficult to achieve with hopping interfer- 
ence paths over third neighbors, even though the square 
root dependence on resistance predicted by hopping the- 
ory is actually closer to experiment than the linear one 
predicted here in Eq. (|45|) . The thermopower and re- 
sistance data by Osinniy et al^ also suggest band edge 
conduction with the usual exponential density of state 
(Eq. {24])) (see Appendix EJ). 



VI. DISCUSSION AND CONCLUSION 

Our analysis suggests that hopping dominates only 
the very low temperature region, and even then, it is 
most likely cluster-to-cluster hopping so that energy 
differences will depend on the cluster magnetism and 
not just on disorder. For temperatures above T ~ lOi^, 
we believe that the Hall coefficient, magnetoresistance 
and thermopower are due to band edge conduction. 
The observed Mott-like laws are most likely caused by 
Coulomb correlations. In band edge transport, any 
value of v between and 1 can be expected. The 
stronger the correlations, the smaller is v. Experiment 
does not give a unique value of v, but it is always 
band edge conduction which is observed. In band 
conduction, the carriers are subjected to incomparably 
stronger magnetoresistance and Hall forces then in 
hopping conduction. This gives rise, in particular, to 
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the anomalous Hall effect. In Appendix [B] we show 
that the conductivity and thermopower data of Osinniy 
et al2& also agree with the model of delocalized band 
transport near the mobility edge. It appears that 
the temperature dependence of the transport data for 
insulating samples depends, to some extent, on sample 
preparation. When the totality of available data in 
the insulating regime of GaMnAs, except those at the 
lowest of temperatures, are examined, it appears to be 
either activated transport or correlation band transport 
which is the dominant mechanism. There are strong 
indications that Mott hopping plays a secondary role 
at this temperature and ferromagnetic concentration 
range. There is, rather, a transfer process by way 
of which a localized hole is annihilated and created 
elsewhere via thermal fluctuations and is a band carrier 
in the intermediate state where it lives long enough to 
be subject to the band Hall forces and spin disorder 
scattering. In other words, the long range band transfer 
mechanism that gives rise to the ferromagnetism in 
the material is also the mechanism that determines the 
transport. The charge transfer process is predominantly 
over the band with a Shklovskii-Efros^i type correlation 
band edge. The picture we propose is similar to that 
of Durst et al— for this regime. It is the dynamics of 
the shared spins, when looked at as real or virtual band 
excitations, which determines the magneto-and Hall 
transport. At very low temperatures there must be a 
region in which hopping at the Fermi level dominates. 
It would be interesting to look for it experimentally. We 
predict that in this range the MR and Hall effect are 
given by Eq. ([25]) together with Eq. {3D]) and Eq. ([15]) 
and are much smaller. Very recently, weak localization 
in GaMnAs has been observe d 53 i 54 . This regime is not 
exactly the one we study here (ej is a somewhat farther 
from the mobility edge in the insulating regime than in 
the weak localization regime) but their results are in 
agreement with our model. 

We have presented a reasonably complete theory of 
transport in magnetic semiconductors of the GaMnAs 
type, when the Fermi level is in the region of localized 
states, in the presence of disorder and magnetism. The 
data available is, unfortunately, sparse and scattered 



and there does not seem to be a complete set of data 
on any one set of samples. The formulae presented 
here have a wide range of applications. The principal 
result of the paper is that we conclude that transport 
in Mn doped GaAs is due to activated transport rather 
than variable range hopping. Finally, the possibility 
of magnetic-cluster-to-cluster hopping and tunneling 
should be taken seriously and investigated in detail. The 
resulting magnetoresistivity and Hall effect properties of 
such systems are not easy to guess. 

Acknowledgments 

The authors gratefully acknowledge financial support 
from the Natural Sciences and Engineering Research 
Council of Canada (NSERC) during this research. P.D. 
also acknowledges support from the Canada Research 
Chair Program. L.-F. A. gratefully acknowledge finan- 
cial support from Pr. Andre-Marie Tremblay during the 
writing of the manuscript. 



APPENDIX A: STEEPEST DESCENT INTEGRAL 

Using a steepest descent approach with Eq. (|3"8j) to 
calculate the conductivity and Hall effect, and neglecting 
the energy dependence, apart from the density of states 
and Fermi function, we have an integral of the type (as- 
sume El = 0) 

7~ / etee" 1 ^ a e > G{e). (Al) 

Je-Ej 

We have changed the sign of the integral from " holes to 
electrons" for convenience. The maximum of the inte- 
grand, neglecting the energy dependence of G is at 

e max ~ [vk B Te v g ] & , (A2) 

provided that e max > \£c — £/|- Clearly this condition 
and ei = are not always satisfied, and depend on the 
material and temperature. When this limit is not valid, 
the integral is controlled by the lower limit and the con- 
ductivity is activated i.e. 



dee kB "' 



G(e) ~k B TG{e f -e. 



,)e (< e /-^>+ e i). 



(A3) 



e-Ef 



When it is valid, the best value can be substituted back, 
and gives a temperature dependence of the form 
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which is a Mott-like law, but without VRH. 



APPENDIX B: THERMOPOWER 

E xx , the energy transported and S, the thermopower 
without interaction effects, are given by Movaghar and 
Schirmacher— : 



Exx 



and 



de 



am 



Pedge{£)Do 



E 



xx, hop 



&E X x 



(Bl) 
(B2) 



In the correlation model with v — 1/3 and for low enough 
temperatures, given by e max > \e c — £/ 



in the insulating regime would be most useful, as they 
should obey Eqs. (|B4[) and (|B5|) and help us under- 
stand the nature of the transport mechanism. Pu et alM- 
have shown that in the metallic regime, the magneto- 
thermopower is anisotropic (depends on the orientation 
of M z to the current). This is due to spin-orbit scat- 
tering. Note that the hopping thermopower, neglecting 
spin flip energy transport, should also depend on the 
magnetic field via the magnetic field dependence of the 
density of states p and the localization length oT x , with 
T = 24 x 2.7a 3 {B z )/[Trk B p{e f (B z ),B z )}. Both spin-up 
and spin-down bands will contribute at the fermi level. 
The localized states have a finite Hubbard U, so the full 
rigorous treatment is not trivial, and needs focused at- 
tention as done in Ref. 58, for the finite U but B-field 
free case. 



APPENDIX C: BAND REGION ABOVE THE 
MOBILITY EDGE 



eT 



(B3) 



This yields a T -1 / 4 law and a magnetic field depen- 
dence, which is part of the Coulomb glass energy e g (B z ) 
(effective Mott To) and not known at present. 



In the exponential band edge model (Eq. I|24p). we 
have 



S = 



s c (B z )-s f (B z ) 
eT 



(B4) 



Finally, in VRH, the thermopower depends on the struc- 
ture of the density of states at the Fermi level. For a 
linear density of states, it is of the for m 28 i 55 i 56 



5' 



k%p(e f )(TT )z 



(B5) 



The form which agrees with the experiments of Osinniy 
et al£2- is the band edge activated law Eq. (|B4|1 . The 
activated form is also what is reported in a-Si:H, where 
the exponential band edge model (Eq. I|24p) applies^. 
The resistance measured by Osinniy et al£* appears to 
also be activated. It seems that we have a material class 
with band conduction, and where the usual amorphous 
exponential density of states model applies. In this ma- 
terial type, Eq. (|A2|) does not apply, and the integral 
Eq. (|A1[) is cut off at fcgT above the mobility edge. The 
transport is now controlled by the mobility edge, and 
Eq. (|B4|) follows. Magnetic thermopower measurements 



The mobility edge contribution to the conductivity, de- 
rived above, is due to the fact that magnetic fields and 
other dephasing effects counteract the quantum Ander- 
son localization process. The B-field moves the mobil- 
ity edge up for holes (down for electrons). But in the 
delocalized region, electrons (holes) are also subject to 
spin-disorder scattering which also causes strong magne- 
toresistance. It may well be that the spin disorder band 
mechanism also contributes when we have the insulating 
scenario, despite the high activation energy. For com- 
pleteness and rigor we therefore need to also examine the 
transport dominated by carriers which are excited below 
(above) the mobility edge in the real band region. When 
the band region starts is not entirely clear but there is 
no reason why spin disorder scattering should not start 
at the mobility edge as well. The band transport should 
therefore include the spin disorder scattering contribu- 
tion as calculated using the Kubo formula^. We know 
that holes excited just below the mobility edge are well 
described by the self-consistent CPA solution 4 .^ so we 
can use it. The contribution from the CPA band con- 
ductivity can be written (e c = s c (B z , M z )) 



e 



de 



dm 

de 



p(e)D C P A (£,B z ), (CI) 



where 



p(e)D CPA (e) = jjJ2j d£kV " [Ihi{Gm(e)}] 2 <5(£-£m)- 

(C2) 
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Also we have 



G k , s (e) 



£ — £k.s — S 



s,CPA 



(C3) 



In Fig. Q] and Fig. [2j we show the calculated curves for 
a parameter range which would correspond to transport 
at the mobility edge. The resistance results shown, mul- 

would give 



tiplied by the activation factor exp 

us a measure of the spin disorder band contribution to 
the magnetoresistance at low temperatures. We note 
that the extended state spin-magnetoresistance is neg- 
ative almost everywhere. This is a result of the fully self- 



consistent single-site approximation which treats both 
the magnetism and the transport on the same footing. 
The intuitive reason is that the material is more ordered 
when the local moments arc aligned. When disordered, 
the holes are scattered from many types of potentials cor- 
responding to the 6 possible S z orientation of the local 
Mn 5/2 moment. When a magnetic field is applied, it 
enhances the order, and this reduces the resistance. Van 
Esch et alZ have shown that the resistance decreases with 
magnetization M Z (B Z , T). We have used their formula to 
analytically simulate the CPA result coming from spin 
disorder scattering in Eq. 
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